Gaussian ellipsoid model for confined polymer systems 
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Polymer systems in slab geometries are studied on the basis of the recently presented Gaussian 
Ellipsoid Model [J. Chem. Phys. 114, 7655 (2001)]. The potential of the confining walls has an 
exponential shape. For homogeneous systems in thermodynamic equilibrium we discuss density, 
orientation and deformation profiles of the polymers close to the walls. For strongly segregated 
mixtures of polymer components A and B equilibrium profiles are studied near a planar interface 
separating A and B rich regions. Spinodal decomposition processes of the mixtures in the presence 
of neutral walls show upon strong confinement an increase of the lateral size of A and B rich domains 
and a slowing down of the demixing kinetics. These findings are in agreement with predictions from 
time dependent Ginzburg— Landau theory. In the case, where one wall periodically favors one of the 
two mixture components over the other, different equilibrium structures emerge and lead to different 
kinetic pathways of spinodal decomposition processes in such systems. 



I. INTRODUCTION 



The tailoring of polymer blends in thin film geometry 
is a problem of active current researchtJ with a wide spec- 
trum of possible technological applications. Of particu- 
lar interest is the design of suitably patterned structures. 
Spontaneous phase separation processes of incompatible 
polymer blends may be used to translate a chemical pat- 
tern on the surface (as e.g. produced by the micro-contact 
printing technique) into a pattern of varying polymer 
compositions.0 Another technique developed recently in- 
volves the alignment jpLpolymer melts in inhomogeneous 
strong electric fields.uTil 

Optimization of these processes requires knowledge 
of the long-time kinetics of polymeric systems on sur- 
faces. Unfortunately, direct numerical investigation of 
this problem by means of molecular dynamics simulations 
is extremely difficult due to the enormous difference be- 
tween typical vibrational times (~ 1 ps) and the charac- 
teristic time scales for diffusional processes of the macro- 
molecules (of the order of seconds). One possible way 
to tackle this problem is to use dynamic self-consistent 
field theoryjj which leads to coupled partial differential 
equations for the monomer densities of the various poly- 
mer components. Another promising approach was sug- 
gested based on the idea that on semi-macroscopic time 
and length scales, polymers canJpe regarded as soft par- 
ticles with an ellipsoidal shapeEl (a different version of 
this idea, using soft spheres instead of ellipsoids, has re- 
cently been proposedlla) . In the work of Ref. 0, a bead 
spring model with repulsive Lennard- Jones interactions 
was used to determine the interaction parameters of the 
soft ellipsoids. Recently, we have set up a similar el- 
lipsoidal model, where the parameters are derived from 
the Gaussian chain model instead]!!] While the underly- 
ing bead spring model is a more realistic basis for de- 
scribing polymers, the Gaussian ellipsoid model (GEM) 
allows one to deal with effective interaction parameters 



that are independent of the chain length, provided the 
self- interaction of the ellipsoids is taken into account. 

In our previous work we have focused on bulk prop- 
erties of polymer melts and mixtures, and verified that 
the model can faithfully reproduce the basic scaling re- 
lations and kinetic properties of polymeric systems on 
large time and length scales. In this publication we ex- 
tend this study to systems in confined geometries. The 
interactions of the ellipsoids with the confining walls are 
modeled by a linear coupling of the monomer densities 
with an external wall potential. The latter is assumed 
to have an exponential form with decay length l w , which 
is numerically very convenient and allows us to tune the 
softness of the wall by changing l w . The detailed pro- 
cedure to account for the wall interactions is presented 
in Sec. 0. Some representative results for homopolymer 
systems are given in Sec. III. In particular, we show den- 



sity profiles as well as wall induced shape and orienta- 
tion effects. The main focus of this work is, however, to 
study phase separation phenomena in binary mixtures 
of two polymers A and B. In Sec. ^ we provide evi- 
dence for surface directed spinodal decomposition waves, 
when the wall prefers one of the components A or B. In 
the case of neutral walls, lateral structure formation is 
observed in accordance with results derived on the ba- 
sis of Ginzburg-Landau type treatments. Moreover, we 
quantify the orientation of the ellipsoids along the phase 
boundaries. 

Finally, in Sec. [v] we consider structured systems, 
where one wall is considered to be chemically patterned 
and the magnitude e w of the wall potential varies along 
one wall direction in a sinusoidal way with periodicity 
iii. This structuring of the surface induces a periodicity 
in the polymer melt. For large film thicknesses periodic 
polymer structures only occur in the vicinity of the pat- 
terned wall. When the film thickness becomes compara- 
ble to L||, the wall induced domain structure propagates 
through the film. 
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II. WALL INTERACTIONS WITHIN THE 
GAUSSIAN ELLIPSOID MODEL 

The Gaussian Ellipsoid Model is described in detail in 
a previous publication.E2l Therefore, we keep the defini- 
tion of the model short and mainly present the necessary 
extensions to describe the interaction of the soft particles 
with the walls. 



A. Free Energy Functional 

The free energy functional of an ensemble of M soft 
particles, corresponding to polymers with iV bonds, is 
given by an intramolecular part, which accounts for the 
possible internal configurations, and an intermolccular 
part, which describes the interaction between the ellip- 
soids. In the presence of boundaries we add a term which 
accounts for the interaction between the particles and the 
surfaces, 
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(1) 



The intramolecular part is determined by the probabil- 
ity -P(S) for a polymer to have the eigenvalues S of the 
radius of gyration tensor S a p (in order S x > S y > S z ), 
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where ks is the Boltzmann constant and T is the tem- 
perature. The intermolccular part is given by 
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where the second term of eq. (|3|) accounts for the self- 
interaction of an ellipsoid. Equation (||) shows that the 
interaction energy is expressed in terms of the monomer 
densities of ellipsoids i and_i (denoted by g • (y) and g'j (y) 
in the laboratory systemllil). High monomer concen- 
trations due to a contraction of a single ellipsoid (self- 
avoidance) or due to a strong overlap of two ellipsoids 
are disfavored energetically (mutual avoidance). The pa- 
rameters e and b 3 are a "contact energy" and a "contact 
volume" , respectively. We denote all reduced energy pa- 
rameters without a hat, e.g. e = e/ksT. The parameter 
b sets the length scale of the model, i.e. 6=1. 

The additional interaction term between the polymers 
and the walls is chosen to be of the form 
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where e w is the strength of the monomer-wall interaction 
and the wall interaction range l w characterizes the hard- 
ness of the wall. The slab thickness is given by L z . Note, 
however, that for the special form of the external poten- 
tial defined in eq. (^J) two parameters e w /ksT = 1 and 
L z are equivalent to two other parameters i w /ksT = e w 
and L z — 2l w In e w , respectively, assuming all other pa- 
rameters to be constant. Without loss of generality we 
therefore set i w /kBT = 1 in the following and define L z 
as our film thickness. This is reasonable, since for l w — > 
the system becomes confined by hard walls at positions 
z = and z = L z . 



An exponential wall potential is frequently used to 
modeL a laterally structureless and purely repulsive 
wallJlj Here, it has the advantage, that the overlap be- 
tween the potential and the monomer densities, that are 
of Gaussian form, are analytically integrable. The for- 
mulas are given in appendix By choosing small values 
for l w one can model hard substrates, while for large l w 
one models soft substrates, that e.g. correspond to poly- 
mer films on polymer brushes &3 We note that the unper- 
turbed monomer densities are used in eq. (^|). Hence, we 
do not a priori build in layering effects which are typically 
observed in simulations of polymer melts close to walls .Eil 
This is in the spirit of the GEM, where the unperturbed 
chain is used as input and all other properties are de- 
termined by interaction terms. In summary, eqs. (0-|^) 
define the thermodynamics of the model for a homoge- 
neous system of ellipsoids in a film geometry. 



To generalize the model to binary mixtures of ellipsoids 
of type A and B in films, we choose an equal interaction 
between polymers of the same type €aa = £bb = e, while 
the interaction between unlike polymers is enhanced by 
a mismatch 6, i.e. eab — e(l + S). The same degree 
of polymerization, i.e. Na = Nb = N, is assumed for 
both components. To account for different interactions 
between the substrates and the two components, we in- 
troduce an additional wall mismatch parameter 5 W , with 
— 1 < 5 W < 1. The interaction parameters e w (l + S w ) and 
e w (l — S w ) refer to the A— and i3-components, respec- 
tively. Note that the mismatch parameters 6 and S w only 
change the strength of the repulsive monomer- monomer 
or monomer-wall interactions. Phase separation is thus 
driven by the disparity of excluded volume interactions, 
and not by attraction. 



The kinetic properties of the GEM are given by a dis- 
crete time Monte-Carlo algorithm, where three different 
moves of the ellipsoids are allowed: Translation of the 
center of mass, free rotation of the principal axis and de- 
formation, i.e. a change of the eigenvalues S,;. For details 
we refer the reader to Rcf. 10. 
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B. Parameters and Geometry in Film Simulations 

The simulations in this work were carried out for el- 
lipsoids corresponding to a chain length of N — 50. The 
number of particles is M — 2000 in the case of homoge- 
neous systems and M = 4000 for binary mixtures. We 
always choose e w = e = 1. 

For binary mixtures a symmetric composition is con- 
sidered, where the fraction of the ^-component is j 'a — 
0.5. For given film thickness L z we determine the lat- 
eral dimension L by requiring a constant "bulk monomer 
concentration" c = M(N+1)/L Z L 2 = 0.85. Tables § and 
U list the geometries of the simulated systems. 

The homogeneous systems are thermalized several 
thousand Monte Carlo steps (MCS) and averaged over 
at least 10000 MCS. In the case of binary mixtures, 
the phase diagrams were determined in the semi-grand 
canonical ensemble. Averages are performed over several 
thousand MCS. 



III. EFFECTS OF HARD WALLS ON 
HOMOGENEOUS SYSTEMS 

In this section we present the results for homogeneous 
ellipsoid systems bound by two walls. As the focus of 
this work will be on binary mixtures in confined geome- 
tries, we defer a more detailed analysis of homogeneous 
systems in a slab to a forthcoming publication. 

As a reference length to compare with, we use the un- 
perturbed average radius of gyration in a bulk system 
with monomer concentration c = 0.85, which is given by, 



cf. Ref. 10 



R G = 3.6 . 



(7) 



reflects the softness of the wall potential, which is also 
responsible for less pronounced layering effects, typically 
observed in polymer melts confined between strongly re- 
pulsive walls.oEj Here, we find the signature of these 
density oscillations only for l w — 0.15, whereas the softer 
wall gives rise to a monomer profile comparable to that 
obtained, for small bulk concentrations in the hard wall 
caseMll 

The monomer concentration profiles for different film 
thicknesses L z are shown in Fig. for the case l w = 0.5. 
For the thicker slabs (L z = 17, 25) a plateau in the pro- 
files is reached at a distance of about 2 Rq from the wall. 
However, the plateau value still varies slightly with the 
system size and is larger than the bulk value c = 0.85, 
that will be obtained for L z — ► oo. For L z < ARc the 
films do not show bulk-like behavior any more. The pro- 
file rather runs through a pronounced maximum c max at 
z = L z /2. 

If L z decreases, starting from L z = 25, c max first in- 
creases, since the reduction of the monomer density close 
the walls becomes more relevant for the overall film be- 
havior, the smaller the thickness is. It has to be compen- 
sated by a higher density in the film center to keep the 
density constant. For L z < 4, however, we observe the 
opposite trend, i.e. a decrease of c max . This is due to the 
fact that the overlap of the wall potentials becomes im- 
portant. This overlap effectively decreases the gradient 
of the wall potential and thus allows the monomers to 
penetrate the wall more deeply. The increased monomer 
density in the wall region entails a reduction of the den- 
sity in the middle of the film in comparison to the thicker 
films. The relative changes of c max are about 20 per- 
cent. In the case of harder walls (l w = 0.15) we found 
a stronger variation and a largest value c max = 1.36 is 
obtained for L = 3. With respect to the binary mix- 



A. Concentration of Monomers and of Centers of 
Mass 



tures studied in Sec. IV we wanted to keep the influence 
of these concentration variations small. We thus chose 
l w = 0.5 in the following. However, very similar results 
could have been obtained for L, = 0.15. 



The results for the concentration of the monomers 

M . 

<z) = Y J Jdxdyg'fy) (8) 

i— 1 

and the concentration of the centers of mass c cm are 
shown in Fig. [I] for a thick film (L z — 25), where bulk-like 
features are found inside the slab. Two wall interaction 
ranges are considered. In the first case, l w — 0.5, mean- 
ing that 2 l w is equal to the average bond length of the 
Gaussian chain molecules, while in the second case the 
interaction range is much smaller, l w = 0.15, modeling 
a harder wall. In both cases, the concentration profiles 
rapidly decay to zero as z — > 0. However, whereas the 
profile of the center of mass, c cm , is zero at z = 0, jus- 
tifying our definition of the film thickness, the monomer 
concentration vanishes only for z < 0. This behavior 



B. Orientation 

In this section we study the orientation of the ellip- 
soids as a function of the distance of the center of mass, 
z, from the wall. 

Let us denote by lZ a f3 the matrix, that defines the ori- 
entation of the principal axis with respect to the labora- 
tory coordinate system. As 1Z a p consists of the normal- 
ized eigenvectors of the principal axis system, J2p T^afj 
1 holds. In an ensemble of randomly oriented ellip- 
soids the average value of all squared matrix elements 
is (T^ad) = 1/3- Thus we expect this value for large dis- 
tances z from the walls. A parallel orientation of one 
of the principal axis a to the laboratory axis j3 corre- 
sponds to {TZ 2 a jj) = 1, while a perpendicular orientation 
corresponds to (T^a) = 0. Hence H^a * s a convenient 
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measure to characterize the orientation. In the following, 
we denote an average of some quantity Q over all parti- 
cles, whose center of mass is at distance z from the left 
wall by (Q(z)). 

Figure [3] shows the results for the orientation of the 
principal axis with respect to the z-axis, i.e. perpen- 
dicular to the walls for different film thicknesses. We 
remind the reader that the system of principal axes 
is chosen such that S x > S y > S z . For thick films 
(L z = 17), we find (Tl 2 az (z)) = 1/3 far from the walls, 
as expected. When approaching the wall the smallest 
principal axis, z, becomes aligned perpendicular to the 
wall, (lZl z (z)) w 1, while the other two principal axes be- 
come aligned parallel £3 This behavior is observed for all 
film thicknesses and is also found in other simulations xB 

The anisotropy in the orientation occurs on a length 
scale 2Rg at most. The range is smaller along the y- 
principal axis than along the ^-principal axis, since due 
to S x > S y the energy cost for a rotation of the y-axis 
in the direction of the z-laboratory axis is smaller than 
that of the ^-principal axis. It is remarkable that the 
anisotropy effect shows only very little dependence on 
the film thickness, even when L z becomes comparable to 
or smaller than Rq. 



C. Deformation 

Figure U shows how the shape of the ellipsoids is altered 
by the presence of the walls. For large distances, we find 
that the averaged eigenvalues (S a (z)} are independent of 
z and coincide with the bulk values. These values, deter- 
mined from the data of Ref. |l^, are S x = 10.1, S y = 2.0 
and S z = 0.71. 

However, the ellipsoids near the walls are distorted. 
The size of the two major axes x and y is enhanced, 
whereas (S z (z)) is reduced. When z increases, (S z (z)} 
approaches the bulk value almost monotonously, while 
(S x (z)) and (S y (z)) run through a minimum before cross- 
ing over to the bulk behavior Jia This means that the ellip- 
soids near the surface are strongly oblate. This influence 
of a repulsive wall on the structure of the chains is also 
confirmed by other simulations (see Ref. [l~5|Jl9| and ref- 
erences therein). As the film thickness L z decreases, we 
find that the maximal values of (S z (z)) become smaller 
for L z < 2R G . 



D. Discussion 

Combining the results of the orientational properties 
and the deformation we find that ellipsoids at a distance 
z > 2Rq are virtually unperturbed by the boundaries 
and exhibit bulk behavior. On approaching the walls 
there is a slight deformation of the two major principal 
axes, accompanied by an alignment of the principal axes. 
This alignment continues monotonously, up to the point, 



where the smallest principal axis of the ellipsoids is ori- 
ented perpendicular to the walls. Near the walls, the el- 
lipsoids are compressed in z-direction and elongated with 
respect to the two major axes, thus exhibiting an almost 
two-dimensional behavior. 

It is interesting to compare these findings with results 
obtained by other authors. Let us mention that the anal- 
ysis of the shape of polymers near walls is usually based 
on "components" of the radius of gyration with respect to 
the laboratory system. These components are the square 
roots of the diagonal elements of the radius of gyration 
tensor, 

RGa = ySota or Rg\\ = ^xx ^VV ' = S zz ■ (9) 

It is clear that 

Sap = TZ ia S 1 TZ 1 p , esp. S aa = TZ^ a S y . (10) 

7 7 

This kind of analysis makes it difficult to distinguish be- 
tween an orientation and a deformation of the polymer 
coils. 

When calculating the quantities (R G ^(z)) and 

(R GJ _(z)) for our model, we find reasonable agreement 
with the results of other simulations (see Ref. |l5|,|l9| and 
references therein), namely an increase of (R G ^(z)} near 

the wall and a strong decrease of (R G± (z)). The results 
for (R G Jz)) are similar to that for (S x (z)) and (S y (z)) 
in Fig. 0, and the results for (R G± (z)) are similar to 
(S z (z)}, except that the extrema near the walls are more 
pronounced. 

Pai-Panandiker et. al@ examined orientational effects 
of the end-to-end vector in a Monte-Carlo Simulation of 
dense polymer films. They found an increasing align- 
ment of the end-to-end vector with reducing L z , but 
averaged their quantities over the whole film. As there is 
a strong correlation between the end-to-end vector and 
the major principal axis x, their results support the pic- 
ture presented above. 

The work of Vliet and BrinkeE] is dealing with a princi- 
pal axis transformation to analyze orientational degrees 
of freedom, for a single self-avoiding walk confined be- 
tween two hard walls. The main conclusion of this work 
is that with decreasing film thickness the principal axis 
of the polymers are aligned due to the confinement, long 
before their shape changes. While their situation is not 
directly comparable to the case of a polymer melt be- 
tween two plates considered here, our results support the 
view that the orientational effect is dominant. This is 
not unexpected, as a reorientation of the polymer only 
constraints on one or two degrees of freedom, while a 
deformation affects the whole chain statistics, therefore 
affecting of the order of N degrees of freedom. 
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IV. DEMIXING PROCESSES IN THIN FILMS 

In this section we consider the effect of two neutral 
(6 W = 0) and planar walls on a binary mixture of ellip- 
soids. Results will be discussed for fixed l w = 1/2. 

A. Dependence of Phase Diagram on Film Thickness 

For temperatures below a critical temperature T c , sep- 
aration into two phases, rich in type A and type B poly- 
mers, occurs. First, we determine how T c depends on 
film thickness L~. Ginzburg-Landau theoryEa predicts 
and simulations^ have shown that 

T c (oo) L z 



T e (co) - T C (L Z ) 1_ 
T c (oo) a LI ' 



L z ^>Nb. 



(12) 



Equations (|ll|,|12j) hold as long as the chain conforma- 
tions are virtually unperturbed, i.e. L z > Rq, and mean 
field theory applies. This requires that the bulk correla- 
tion length £ (characterizing the interfaciajL-svidth close 
to T c ) is smaller than the crossover lcngthEa £ x cx Nb, 
which determines the change from Ising-type to mean 
field behavior. While for L z <C Nb the surface interac- 
tion is dominant [eq. (jll|)1, it is the finite size of the film, 
which determines the shift of T c for L z 3> Nb [eq. (Q)]. 
Since N — 50 and 2 < L z < 25 in our case, we expect 
eq. ( pi] ) to apply at least for L z > Rq — 3.6 [eq. (^)] . 

For the sake of efficiency the simulations are carried 
out in a semi-grand canonical ensemble£j i.e. the total 
number of particles is fixed, while the fractions Ja arL d 
fs of A and B polymers can fluctuate. Figure ^ shows 
the order parameter (\/a ~ /b|) as a function of the mis- 
match parameter 6 for various L z . 

In order to test the expectation that our model obeys 
eq. (p"l|), we simply estimated 6 C oc 1/T C from an extrap- 
olation of the steepest descent of the coexistence curve to 
zero. No detailed finite size analysis of the lateral dimen- 
sion L was attempted. Such an analysis does not seem 
to be really necessary, since T C (L,L Z ) monotonously ap- 
proaches the asymptotic limit T C {L Z ) for large enough 
Lx3 For the L- values used in this study (see Tabic fil), 
inspection of Fig. 5 of Ref. ^ suggests that qualitatively 
reliable results should be obtained from our estimates of 
T c . In fact, we see that T c decreases with L z . The in- 
set of Fig. U shows that T C (L Z ) — T c (oo) cx —1/L Z , in 
agreement with eq. (|ll|). 



B. Orientation at Phase Boundaries 

To study the influence of phase boundaries on the in- 
ternal degrees of freedom of the ellipsoids, we prepared a 
fully phase separated system with periodic boundary con- 
ditions in all directions L = L z = 62.14. The A-ellipsoids 



are randomly distributed in the half-space z = (0, L z /2) 
and the -B-ellipsoids in the half-space z = (L z /2, L z ). 
We chose S — 0.2, which corresponds to the strong segre- 
gation regime (cf. Fig. |^) for thick films. The system was 
thermalized for 2000 MCS, a time much larger than the 
autocorrelation time of the internal degrees of freedom 
(shape and orientation) . Moreover, this time is sufficient 
for an ellipsoid to diffuse typically over the distance from 
one interface to the center of one of the two phases. All 
averages were recorded over additional 18000 MCS. 

In Fig. ^a the different center of mass densities for 
both components Cc m> A, c cm , s , and c cm = c cm:A + c cm ,s, 
normalized by the chain length (N = 50) are displayed. 
We see that the total density c cm of the ellipsoids be- 
comes smaller at the interfaces in order to reduce the 
interaction between the mixture components. An addi- 
tional effect occurs in Fig. ^b, where the squared matrix 
elements {Tt 2 az ) are show n. Comparable to the behavior 



near a wall (cf. Sec. [II B) the ellipsoids align their x-axis 
(largest principal axis) parallel to the interface, while the 
z-axis (smallest principal axis) tends to be oriented per- 
pendicular to the interface. However, this effect is much 
less pronounced than for walls. In particular, there is 
almost no effect on the orientation of the y-axis. 

Finally, Fig. ||c characterizes the deformation of the 
ellipsoids via the eigenvalues S a that are normalized by 
their bulk values. Near the interfaces the ellipsoids are 
compressed with respect to all axes. The relative com- 
pression increases with the length of the principal axis, 
meaning that the asphericity becomes reduced. This ef- 
fect was also observed in extensive Monte Qarlo sim- 
ulations of the bond-fluctuation lattice modelEil and is 
analogous to the reduction of the radius of gyration in 
a strongly segregated mixture, where the polymers of 
the minority component shrink to reduce energetically 
unfavorable contacts with the the surrounding majority 
phased 



C. Spinodal Decomposition Between Neutral Walls 



Studies based on time-dependent Ginzburg-Landau 
theoryEH have shown that for spinodal decomposition in 
thin neutral films lateral decomposition waves of a char- 
acteristic wavelength fen m characterize the domain pat- 
tern at short times, where the concentration inside the A 
and B rich domains has not yet reached the equilibrium 
value. These studies predict that feii m becomes smaller 
with decreasing L z , when L z < A m , where X m denotes 
the bulk demixing length. Furthermore, if L z < A m , the 
growth rates of the decomposition modes become smaller 
so that the overall demixing process slows down. These 
predictions will be tested here for the GEM. 

We consider a system confined between two planar and 
neutral walls (6 W = 0) that is prepared at infinite tem- 
perature at a time t < 0. All orientations and particle 
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positions are random and the distribution of the eigen- 
values is given by those of free Gaussian chains. At t = 
we quench the system into the spinodal region of the 
phase diagram by setting e = t/ksT — 1 and 6 = 0.2 
(for L z = 2 we used 6 — 0.25). The dynamics of the 
system is specified by Monte-Carlo moves of the indi- 
vidual ellipsoids corresponding to translations, reorien- 
tations and deformations. The detailed procedure is de- 
scribed in Ref. [l0| 

To quantify the demixing process of the films, we de- 
termine the lateral intermediate scattering function 



length and the slowing down of the demixing kinetics for 
L z < A m . This qualitative agreement with the results 



7|| (fe|| ,z,t) = (p A \\ (k|| , z, t)p A \\ (-k|| , z, t)) 



c 2 A S(k v 



(13) 



where p A \\ (k|| , z, t) is the Fourier transform of the to- 
tal monomer density g' A (y,t) = X)i=i Qi A^' *) °^ * ne 
A-chains with respect to the lateral coordinates x, y. 

1 /2 

The power-spectrum depends on fey = (k 2 + ky) only, 
therefoija.it is possible to perform an additional circular 
averageta in the k x , fc y -space. An explicit expression for 
Pa\\ (ky , z, t) is given in appendix [b]. 

Figures 0a-c show I»(k\\,Z,t) for a slab of thickness 
L z = 6 and various fixed times. As one can see from 
these figures, Iu(k\\,z, t) exhibits a maximum that cor- 
responds to a characteristic size 2ir/k^ m of the A and 
B rich domains in the lateral direction. With increasing 
time the maximum grows and shifts to smaller values of 
fc|| , reflecting lateral domain coarsening. 

To quantify the speed of this domain coarsening and to 
examine its dependence on film thickness L z , waxonsider 
the first moment of hi (fc|| , z, t) averaged over zjij 



*i,||(*) 



1 



n / °° dk\\ I\\ (fc|| , z, t) 



(14) 



This quantity is a measure of the average lateral domain 
size as a function of time. The results are given in Fig. || 
for different values of L z . 

In our system the bulk demixing length A m can be es- 
timated by X m w 27r/fci(t Q ) « 14, cf. Fig. 15 of Ref. [h] 
(to = 500 MCS). From this we expect that the effects 
predicted by time-dependent Ginzburg-Landau theory 
should occur for L z < 14. 

Indeed, we find that (*o) determined at t — 1000 
MCS shrinks with decreasing L z , cf. Fig. B. Moreover, 
in the thinner films the coarsening sets in later, which 
is illustrated by the fact that curves at equal times have 
a steeper slope for thick films than for thin ones. This 
supports the result that the growth rates of the decom- 
position modes decrease. We note that these effects are 
not due to a shift of the critical temperature T c with 
L z . We carried out additional simulations with a mis- 
match 6 — 0.3 for L z = 2, 3, 4 and found that the lateral 
demixing length depends only weakly on the mismatch 
parameter at fixed film width. 

To summarize this section, our results support the pre- 
dictions of Ref. Eft i.e. the increase of the lateral demixing 



of Ref. 27 is obtained without any special wall potential 
or tuning of the parameters, supporting the view that 
the effects occur quite generally. However, a suppression 
of spinodal decomposition below a critical film thickness 
L c , as predicted in Ref. |?], is not found. This is due to 
the fact that noise terms in the kinetic equations of the 
theory have been neglected. 



D. Surface Directed Spinodal Decomposition 

In the case, where one wall preferentially attracts one 
of the two components, surface directed spinodal decom- 
position wavesta (SDW) form when a homogeneous mix- 
ture is quenched into a thermodynamically unstable part 
of the phase diagram. However, in thin films the SDW 
might be "suppressed" by the confinement £3 

Here, we investigate the case, where both walls attract 
ellipsoids of component B. Figure ^| shows the emerging 
monomer concentrations cg(z), averaged over the lateral 
dimensions. For rather large slabs (L z = 50) symmet- 
ric SDW's form, cf. Fig. |^a. Their wavelength at short 
times is given by the bulk demixing length X rn w 14, 
and at large times increases due to domain coarsening. 
Finally, at t = 100000 MCS a B-A-B stripped pattern 
has formed in the z-direction. In the case of thin films 
(L z = 12.5, cf. Fig. ^|b), this kind of pattern emerges 
already for short times (t = 1000 MCS). For L z < \ m 
the propagation of the SDW is suppressed, and only a 
fast segregation of the B-components at the walls occur. 
For even thinner films (L z = 6) the two maxima seen in 
Fig. |^b, are less pronounced and merge into one in the 
case of L z = 3. 



V. STRUCTURED SURFACES 

In this section we study polymer mixtures between two 
walls where one of the two surfaces is chemically pat- 
terned. One part of the patterned wall attracts compo- 
nent A, while the other part attracts component B. The 
second homogeneous wall is neutral with respect to the 
two components. The surface pattern is a periodic ar- 
rangement of equally spaced stripes, cf. Fig. [l(], modeled 
by a surface potential 



V(z) = e w exp ( - — 

Z In) 



1± 5 W cos 



2iry 



(15) 



The positive sign refers to component A, the negative 
to component B. As for homogeneous walls, the over- 
lap integral [eq. (||)] can be calculated analytically, see 
App. |X| In the following we consider the case l w = 0.5 
and S w = 0.5, where both components are strongly fa- 
vored by the respective parts of the wall. Our focus will 
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be on the influence of the film thickness L z and of the 
pattern periodicity Lu on the demixing process. 



A. Structural Phase Diagram 

From a technological point of view it is important to 
explore regions in L z — Lu space, where the pattern of the 
surface can be translated into a corresponding polymer 
structure that pervades the whole filmBEl 

We carried out simulations in a semi-grand canonical 
ensemble for different combinations of L z and Lu sum- 



marized in Tab. III. As we require L y to be an integer 
multiple of Lu and work with a fixed number of particles 
M = 4000, we adjusted L x so that the monomer den- 
sity (or volume) is constant. The quantities of interest 
are determined after an equilibrated stationary domain 
pattern has evolved. 

To characterize the equilibrium patterns we recorded 
the center of mass density c cm . Two representative re- 
sults for a slab of thickness L z = 6 are given in Fig. It 
shows the normalized center of mass density (JV+1) c cm 
as a function of one laboratory axis (y or z), while the 
density profile was integrated over the respective remain- 



ing axes. For Lu = 25, Fig. 11a, the A and B ellipsoids 



are separated in a periodic manner due to the structure 
of the wall. As one can see from the inset, this structure 
propagates through the whole film ( "full structure" ) . By 
contrast, for Lu = 6, Fig. pl]b, the periodic behavior is 
confined to a small region near the structured surface 
(cf. the inset). In this case, small cylindrical caps of 
one phase are present at the favored parts of the surface 
("partial structure"). These situations, full and partial 
structures, are schematically depicted in Fig. 

A systematic study of the emergingstructures as a 
function of L z and Lu is shown in Fig. 12. The (Lu,L z )- 
plane is divided into two regions, where systems with 
small Lii and large L z are only partially structured, while 
systems with large Lu and small L z are fully structured. 
To explain this behavior, we associate with the half pe- 
riods of the cosine term in eq. (|15[ ) two homogeneous 
surface regions 7 and <5, cf. Fig. [10[ . The contact angles 
in these regions are 



COS0, 



■ cos 85 = . (16) 

OAB &AB 



Near the j-S interface the contact angle can assume 
any value between the two limiting cases 7 and @<5.E2l 
Requiring without loss of generality that ca < cb, the 
contact angle is restricted toO<G<7r/2. A simple 
consideration of the area of the interfaces in both cases 
(a) and (b) given in Fig. [l(] yields a difference of the free 
energy 



AF = cjabL x L v 2L 



2 sin9 11 



(17) 



The transition between both structures occurs, if AF 
0, resulting in 



hi = 6 
L|| 4 sin 6 



(18) 



The crossover regime between the two structures in the 
L z —Lu diagram is between the two lines with slope 1/4 
(for 6^0) and slope tt/8 (for 9 -> tt/2). Both lines 
are marked in Fig. |l2| for comparison, showing agreement 
with the simulations. 

This result is also confirmed tw experiments of an A/C 
demixing process in thin filmsHJ In the experiment, the 
patterned substrate was realized by physisorbing the B- 
component of an ABC triblock copolymer microphasc 
separated brush. This leads to a checkerboard like A/C 
surface structure, instead of the striped structure consid- 
ered above. However, the different geometry yields only 
a different constant in the righthand side of eq. (|l^). 



B. Pattern Directed Spinodal Decomposition 

In this section we consider the kinetics of a system sim- 
ilar to the one described in the previous section, but for 
fixed fraction /a — 0.5. The initial conditions are as in 
Sec.lfVC. 



The intermediate scattering function 



Iy(k y , z, t) 



dk x J||(k||, z, t) 



(19) 



with lu (k|| , z, t) given in eq. ( |l3| ) characterizes the domain 
formation with respect to the pattern-direction. 

Figure [l^ shows I y (k y , z, t) for a thin slab of thickness 
L z = 6 and pattern periodicity L|| = 25, corresponding 
to a full structure in equilibrium, at two times. At early 
times (f = 1000 MCS, see Fig. 13a,), a dominant structure 
arises near the patterned substrate on the length scale 
ky l = L||/27r. This structure propagates in z-direction 
until it reaches the opposite wall (t = 10000 MCS, see 
Fig. |l3|b). At later times, this structure is stable, as one 
would expect from Fig. [12| We determined no further 
change in Uk y ,z, t) until t = 30000 MCS. 

Figure [14] shows I y (k y ,z,t) for a film of thickness 
L z = 25 and pattern periodicity Lu = 12.5, correspond- 
ing to a partial structure in equilibrium. In this sys- 
tem the pattern-induced demixing is much faster than 
the lateral spinodal decomposition. The pattern-induced 
structure for small values of z leads to a much larger cor- 
responding peak in I y (k y , Z,t), Fig. |l4|a, than the disor- 
dered lateral domain pattern emerging for larger z. The 
ordered periodic structure propagates into the film over 
a finite distance only [z w 6, see Fig. |l4|b). For later 
times (t = 30000 MCS), the amplitude of the lateral de- 
composition waves is comparable to the amplitude of the 
pattern-induced waves. However, the peak in I v {k y , z,t) 
occurs at smaller values of k y . Thus, the periodic domain 
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structure is stable near the patterned wall, but further 
away the domain pattern coarsens. 

Such a behavior, i.e. a rapid formation of a patterned 
equilibrium structure in thin films, but a lateral do- 
main coarsening afterwards in thick films, has also been 
observed in numerical treatments of the Cahnp 
equation with appropriate boundary conditions^ 
in experiments. 

In both cases considered so far the periodicity Lu was 
comparable to or larger than the bulk demixing length 
A m . The question arises of what happens for Lu <C A m . 
In this case, we expect that the patterned surface induces 
structures, which are too small to drive spinodal decom- 
position because the interfacial energy of such structures 
is higher than the free energy gained by decomposition. 

To investigate the interplay between the two length 
scales L|| and A m _we consider a slab of thickness L z — 
12.5 and vary iyJ^J Figure [l5] shows corresponding lat- 
eral structure factors I y (k y ,z,t) for t = 1000 MCS and 
for Lu = 6, 12.5, and 25. In the case of Ln = 6 <C A m , 
Fig. |l5|a, the peaks corresponding to the surface in- 
duces structure and the lateral modes are comparable. 
The demixing process is not dominated by the surface- 
induced pattern. Also, the induced structure develops 
much more slowly compared to the case of Lu — 12.5 ~ 



A m , shown in Fig. |15|b. In this case, the peak corre- 
sponding to the ordered structure is a factor of about 10 
higher, and the lateral demixing modes are thus almost 
invisible. For Lu = 25 ^> A m a double peak structure 
emerges. This can be interpreted in terms of a pattern 
directed spinodal decomposition wave, i.e. if one looked 
at the concentration fluctuations ScA(y) = f?>i(y) — C A m 
the yz- plane and labeled positive and negative regions 
of ScA{y) black and white- respectively, one would see a 
checkerboard-like patternO 

To see more directly how in the last case such a pattern 
propagates into the film, we plotted in Fig. |l^ the inter- 
mediate scattering function 7j,(27r/L||, z, t) as a function 
of z for various times t for a system with L z — 25 and 
Lu = 50. This picture describes how the surface-induced 
periodic domain patterns grow in perpendicular direc- 
tion. Since our film thickness is limited, we did not try 
to quantify the propagation of this kind of kinetics fur- 
ther. For very large film thicknesses we would expect the 
usual Lifshitz-Slyozov domain growth to occur at large 
times. 



VI. SUMMARY 

The Gaussian Ellipsoid Model0 (GEM) is based on 
the fact that the most probabJe^ihape of a polymer is 
not spherical, but ellipsoidal.EjEj This ellipsoid resem- 
bles a flattened American football, the eigenvalues S of 
the three principal axes being all different. In the GEM, 
the eigenvalues and their distribution are calculated for 
a Gaussian chain. Thus, the model should faithfully rep- 



resent the large scale chain structure in a 8-solvent or in 
a melt. In order to extend the model to dilute and semi- 
dilute solutions in good solvents excluded volume interac- 
tions have to be added. This is achieved by determining 
the local monomer density in the Gaussian ellipsoid. A 
strong overlap of the density in the ellipsoid or between 
different ellipsoids entails a large energy penalty which 
mimics the self- and mutual excluded volume interaction 
of soft polymer coils. 

In the bulk, the model reproduces well established 
properties on large length and time scales. lj For ho- 
mopolymers, it yields the correct crossover scaling from 
dilute to semi-dilute solutions and the correlation hole 
in the melt. For symmetric homopolymer mixtures, one 
(approximately) finds a linear dependence of the critical 
temperature on chain length, a contraction of the chain 
size of the minority component, and a Lifshitz-Slyozov 
behavior for the late stages of spinodal decomposition. 

These results prompt us to test the model if interfaces 
are present. These interfaces are either self- generated, as 
it is the case in the strong-segregation regime of binary 
polymer mixtures, or imposed externally by the intro- 
duction of two repulsive walls on opposite sides of the 
simulation box. The main results of this extension may 
be summarized as follows: 

(1) In the vicinity of a repulsive wall the ellipsoids can 
no longer adopt all possible orientations. They become 
aligned and distorted. The two longer axes are stretched 
and oriented parallel to the wall, whereas the short axis 
shrinks and orients perpendicular to it. These wall- 
induced perturbations of the chain's structure cross over 
to the isotropic bulk behavior for large film thicknesses 
(i.e., L z > ARq)- These results parallel those obtained 
from less coarse-grained simulation models (Refs. p^JT^ 
and references therein). 

(2) Chains close to the interface in a strongly segre- 
gated polymer mixture behave similarly. In order to 
reduce unfavorable interactions the largest axis of the 
ellipsoids orients parallel to the interface and shrinks 
in size so as to increase the number of favorable self- 
contacts inside the chain. These results qualitatively 
agree with those obtained in a detailed Monte Carlo 
study of the phase behavior of polymer blends by the 
bond-fluctuation model.Ej 

(3) If a binary mixture of ellipsoids between two neu- 
tral walls undergoes spinodal decomposition, a lateral 
demixing length An TO is found, which increases for films 
thinner than the bulk demixing length (L z < A m ). These 
results are in affpeement with time-dependent Ginzburg- 
Landau theory.cH Experimental studies of this effect have, 
to the best of our knowledge, not been performed yet. 
However, both the present simulations and theory lead 
to similar results and suggest that such studies should be 
promising. 

(4) In the case where one of the two walls is structured, 
two different equilibrium patterns are found. Either the 
structure propagates through the whole film, leading to 
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a striped demixing pattern, or the surface-induced struc- 
ture is confined to the vicinity of the wall. Both types 
of pattern are separated by a straight line in the L z -L\\ 
space, as can be explained by a consideration of surface 
energies only. Different equilibrium patterns lead to dif- 
ferent kinetic pathways for the spinodal decomposition in 
such a film. In the case of an equilibrium pattern with full 
structure, pattern directed spinodal decomposition is the 
dominant process. The periodic ordered pattern freezes 
after reaching the opposite wall. On the other hand, the 
penetration depth of the pattern directed spinodal wave 
remains finite and lateral domain coarsening takes place 
in the late stages of decomposition. In the case of small 
pattern periodicity (L\\ < A m ), pattern directed spinodal 
decomposition is less pronounced, as the pattern induces 
unfavorable structures with too large interfaces. 

In summary, the presented results suggest that the 
GEM is a promising model to study large length and time 
scales not only in the bulk, but also in the presence of in- 
terfaces. To test the model further we envisage to analyze 
the dependence of the structure of a .thin polymer film 
on chain length. Several simulationsOEj have already 
been devoted to this problem, with which we can com-, 
pare our results. Furthermore, recent experimentscda 
investigated the dependence of the radius of gyration on 
film thickness and obtained diverging results. We hope to 
contribute to this current discussion by the simulations 
envisaged. 
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APPENDIX A: CALCULATION OF 
ELLIPSOID- WALL INTERACTION 



x(y) = 



■ exp 



4e 



Ca(y)* 



c Q (y) = E H a p(yp ~wp), a a ^R a a a and 



tp(y 3 ) = exp ( ±-y- 



(A2) 



We adopted the notation and definitions of Ref. [lC]. An 
analogous calculation yields the same result for the inte- 
gral as in eqs. (C12-C14) of Ref. [l^, except for the values 
of A a p,B a and C, that in the present case read 



•Act/3 — £*a/3 j 

C = ^ C a pW a Wp . 

a, (3 



(A3) 
(A4) 

(A5) 



The ne gative sign in front of the Kronecker-delta in 
eq. (A4) corresponds to the positive sign in eq. ( |A2| ) and 
vice versa. 

In the case of a patterned substrate, we define lu = 
L||/47r and write for the first part of the wall potential 



V(y 3 ) = e w exp 



V3 

21 



1 ± <L, cosh 



m 

21 



(A6) 



Thus the integration of the respective exponential func- 
tions can be performed as above. Only the form of B a 
changes, which now reads 



B n 



T — 2 CfiaWp , 





(A7) 



where the negative sign ref ers to the first imaginary ex- 
ponential function in eq. ( |A6| ), while the positive sign 
refers to the second one, respectively. 

We note, that in our implementation the interaction 
between the walls and the ellipsoids is never cut off, in 
contrast to the interactions between the ellipsoids .E2l 



The integrals in eq. (||) with a potential of the form 
eq. (^|) and monomer densities, that are essentially sums 
and products of Gaussian functions, cf. eqs. (25a-c) in 
Ref. [h], can be cast into the more general form already 
treated in the appendix of Ref. [l(J. The six resulting 
integrals of eq. (|5|) are of the form 



APPENDIX B: CALCULATION OF THE 
LATERAL FOURIER TRANSFORM OF THE 
DENSITY 



The calculation of 



+oo 



d y x(y)^(ya) ■ 



where the functions x(y) an d "0(2/3) are 



(Al) 



pw ( k n . y 3 ) = E j_J 2 y\\ e " lk|! ' y " e'M > ^) 

(Bl) 







is analogous to the derivation of p(k) given in appendix 
D of Ref. nG. Therefore we only present the result, 



P||(k||,?/3) = 



1 



1 



1 M 

-T 



1 



2tt (2 + ci)ctict 2 CT3 M j-^ J s (i) s ^) s ^ 



(0 



exp 



V "=1,3=1 / 



exp 



exp 



2 sir 



/9=1 



2 cos 



i8=l 



wr 1 



x < exp 



kaKiaUiySl 1 



with 

y a/3 



?7 /v 7«' v 7/3 



G 



(0 
3/3 



7 

033 +Saa and 



The quantities 72.^1 and rJ, IJ denote the rotation ma- 
trix and the center of mass of the «'th particle respec- 
tively. The constants a± , 02 , 03 , c \ and u\ are given in 
Tab. V of Ref. [u]. Using eqs. (p2}p5|) it is straight- 
forward to calculate the lateral intermediate scattering 
function, eq. (13). 
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TABLE I. System sizes for homogeneous films. 
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TABLE II. System sizes for films of binary mixtures. 
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TABLE III. System sizes for films of binary mixtures on 
structured surfaces. 
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FIG. 1. Monomer density c(z) (upper curves) and center 
of mass density c cm (z) (lower curves) for wall potentials with 
l w — 0.5 (open symbols) and l w = 0.15 (full symbols). The 
film width is L z = 25. 
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FIG. 2. Monomer density as a function of distance z from 
the wall for l w = 0.5 for various film thicknesses L z . 
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FIG. 3. Lateral averaged squared matrix elements 
(^^(z)) are shown characterizing the orientation of the ellip- 
soids with respect to the z-axis of the system. The alignment 
of the principal axis corresponding to the largest eigenvalue 
is given in (a), to the second eigenvalue in (b), and to the 
smallest one in (c). The expected value 1/3 for totally ran- 
dom oriented ellipsoids is marked in as a straight line each 
time. The results are shown for different film widths L z and 
l w = 0.5. 
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FIG. 5. Coexistence curves for different film thicknesses 
L z = 25(D), 12.5(B), 6(o),4(»), 3(A), 2(A). The extrapolated 
critical temperatures as a function of 1/L Z are given in the 
inset. 
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FIG. 4. Lateral averaged eigenvalues (S a (z)) as a function 
of the separation z from the wall. The expected bulknval- 
ues from simulations with periodic boundary conditionaij are 
shown additionally. The results are given for different film 
widths L z and l w — 0.5. 
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FIG. 6. Normalized center of mass density (a), squared ro- 
tational matrix elements (b) and eigenvalues (c) as a function 
of distance z from the interfaces for a separated A-B mixture. 



FIG. 7. Lateral intermediate scattering function lu (ku , z, t) 
for a film of thickness L z = 6 shown at various times (a) 
t = 3000 MCS, (b) t = 10000 MCS and (c) t = 50000 MCS. 
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FIG. 8. Film averaged first moment of the lateral inter- 
mediate scattering function ki t ^(t), eq. ([l4|), as a function of 
time for different film thicknesses L z . 
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FIG. 10. Two possible structures on a patterned substrate. 
Upper figure (a): The polymer film is only patterned near the 
surface ("partial structure"). Lower figure (b): The polymer 
structure induced by the surface penetrates the film ("full 
structure"). 
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FIG. 9. Laterally averaged monomer concentration cb(z) 
of the i? -ellipsoids, for different times, in a film of thickness 
(a) L z = 50 and (b) L z — 12.5. The mismatch wall interac- 
tion is 8 W — 0.5 for both walls. 



FIG. 11. Results for the normalized center of mass density 
Can as a function of y and z (shown in the respective inset) 
in the case of (a) full, and (b) partial structure. The film 
thickness is L z — 6 in both cases and the periodicity of the 
structure is (a) L\\ = 25 and (b) L\\ = 6. 
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FIG. 12. Regions in the L z -L\\ plane, where the domain 
pattern in the slab is fully structured (■) and only partially 
structured (□), in double logarithmic representation. The di- 
viding lines are given by the two limiting cases of eq. (CLS[) , 
L z = 1/4 Z/ 1| (straight line) and L z = 7r/8L| (dashed linejT^ 
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FIG. 14. Intermediate scattering function I y (k y ,z,t) for a 
system of L,| = 12.5 and L z = 25 at times (a) t = 1000 MCS 
and (b) t = 30000 MCS. Lateral domain coarsening occurs 
away from the vicinity of the patterned wall. Near the wall 
the structure changes only weakly. 



FIG. 13. Intermediate scattering function I y (k y ,z,t) for a 
system of L\\ = 25 and L z = 6 at times (a) t = 1000 MCS 
and (b) t = 10000 MCS. The patterned structure is stable 
throughout the whole film. 
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FIG. 15. Intermediate scattering function I y (k y ,z,t) for 
early demixing times t — 1000 MCS and film thickness 
L z = 12.5. The pattern periodicity is (a) Ly = 6, (b) 
L|| = 12.5 and (c) L\\ = 25. 



18 



